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Abstract 

We study the tunneling of a small ensemble of strongly repulsive bosons in a one-dimensional 
triple-well potential. The usual treatment within the single-band approximation suggests suppres- 
sion of tunneling in the strong interaction regime. However, we show that several windows of 
enhanced tunneling are opened in this regime. This enhanced tunneling results from higher band 
contributions, and has the character of interband tunneling. It can give rise to various tunnel- 
ing processes, such as single-boson tunneling and two-boson correlated tunneling of the ensemble 
of bosons, and is robust against deformations of the triple well potential. We introduce a basis 
of generalized number states including all contributing bands to explain the interband tunneling, 
and demonstrate various processes of interband tunneling and its robustness by numerically exact 
calculation. 
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I. INTRODUCTION 



Ultracold atoms in optical lattices represent a rapidly-growing field of research with 
applications such as quantum simulators j4| of strongly correlated systems, quantum phase 
tra.s,ti„„s. such as the superfluid-Mo|t msulator transition ffl. or quantum computing 
and quantum information processing |8|, |9[. All these applications are rooted in the possibil- 
ity of an extensive control of the underlying setup: the geometry and strength of trapping 



potentials can be tuned by the laser beams forming the lattice potential 
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I0I-I12I, while the 



15] . Particularly in quasi 



interaction strength can be changed via Feshbach resonances 
one- dimensional systems, where the transversal degrees of freedom are energetically frozen 
by a tight transversal confinement, the effective interaction strength can be tuned by the 
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2Q|. As a 



confinement strength, leading to so-called confinement induced resonances 
hallmark example, this has led to the realization of the Tonks- Girardeau gas, a highly corre- 
lated one-dimensional bosonic state where the strong repulsive interaction induces fermionic 



properties into bosons 
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2^. 



The double well is the simplest case of a multiwell potential, and a plethora of phe- 
nomena have been observed and analyzed in great detail for this system. Loaded with a 



Bose- Einstein Condensate (BEC), the doub 



the superconducting Josephson junction 
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e well manifests itself as a bosonic analogue of 



26] . In the low-interaction regime, the coher- 



ent phase coupling of the BEC in neighboring wells dominates the tunneling properties, and 
gives rise to population transfer such as Rabi oscillations and 7r-mode oscillations 27|, while 



in the somewhat stronger-interaction regime, the nonlinearity introduced by the interaction 
dominates, and the BEC can become self-trapped 
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29] . On the other hand, for few-body 



systems, the microscopic counterpart of BECs, the tunneling between two neighboring wells 



also takes place in a corre. 
interaction regime 
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ated manner, exhibiting pair tunneling in the low- to intermediate- 
3l| . and turning to self-trapping for stronger interactions 31|, which 
indicates the intrinsic correlation between micro- and macro-ensembles of ultracold atoms. 
The generalization from double well to multiwell systems can provide a bottom-up under- 
standing of mechanisms operating in the infinite optical lattice. For instance, the study of 
the ground state properties by increasing the interaction strength in the double well and 
multiwells with the well number up to seven shows a common behavior of vanishing inter- 



well correlations 



which is related to the superfluid-Mott insulator transition in the 



infinite optical lattice. 

A straightforward natural extension of the double well is the one-dimensional triple well. 
In the triple well, ultracold atoms also show correlated tunneling and self-trapping in the 
lower- and stronger interaction regime, respectively 35|, and present at the same time novel 
properties such as stationary tunneling, in which the loss of coherence between bosons leads 
to a steady state 36|. Reference 37| provides an extensive study of the long-time evolution 



of the dynamical properties of micro- and macro-ensembles of bosons in the triple well. All 
these studies may generalize the phenomena in the double well to the optical lattices, and 
bridge the gap between the micro- and macro-systems, as we can see the triple well as a 
protype of optical lattices. Moreover, similarly to the analogue of the double well to the 
Josephson junction, the triple well is the minimal system which can model the source-gate- 
drain junction, and draws lots of attention from the perspective of atomtronics. A variety of 
proposals to achieve controllable atom transport based on the triple well have been presented 
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A widely used approximation in the study of cold atoms is the single-band approxima- 
tion, which ignores the higher band contribution. This single-band approximation, working 
successfully in the low interaction regime, finds difficulties in the strong interaction regime, 
where the strong correlation between bosons suppresses the coherence, and consequently 
leads to phenomena such as the fragmentation of bosonic ensembles (4l|, fragmented pair 
Tunneling and eventually breaking up of the tunneling pairs in the To nks- Girardeau regime 
42|. Efforts have been done to extend the studies to the strong- interaction regime, by 
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explicitly taking into account the higher band contributions. It has been shown that higher 
band effects can give rise to effective many-bod y in teractions 



Bloch oscillation 



45| as well as Mott transition 
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46-48 
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44| . and can modify the 



In this work, we study the dynamical properties of a few-boson system in a one- 
dimensional triple-well with the numerically exact Mult i- Configuration Time-Dependent 
Hartree Method (MCTDH) 49N5l|. In the strong interaction regime, where it is common 
sense that tunneling in the multiwell is suppressed, we observe that multiple windows of 
enhanced tunneling are opened on top of the suppressed-tunneling background. These tun- 
neling revivals result from the resonant coupling of the initial state to particular higher 
band states, and this enhanced tunneling is interband tunneling. Various tunneling patterns 
such as single-boson and correlated tunneling can be realized in corresponding windows of 



enhanced tunneling, which could be used to transport cold atoms in a controllable way. 

The present work is organized as follows: in Sec. II, our setup and the MCTDH method 
are introduced. In Sec. Ill we discuss the construction of generalized number states which 
include interactions and are ideally suited for the analysis of the interband tunneling. In 
Sec. IV we present our results and their discussion and interpretation, and demonstrate in 
particular the interband process. Sec. V contains the summary and conclusions. 



II. SETUP AND COMPUTATIONAL METHOD 



A. Setup 

In this work, we explore the correlated quantum dynamics of a few-boson system confined 
to a one- dimensional triple well. The Hamiltonian takes on the appearance: 

N ^2 ^ 1 

i=i i=i 

Here Vtr{x) = Vosm^(/tx) models the ID triple-well confinement, where for convenience we 
apply hard- wall boundary conditions at x = ±37r/2K (k is the wave vector of the laser beams 
forming the optical lattice) so as to confine the bosons to three adjacent wells, as schemati- 
cally shown in figure 1(a). Experimentally the triple well potential can be realized, e.g., by 
a bichromatic optical lattice, formed by laser beams of different frequency, or a harmonic 
trap superimposing an optical lattice. The shape of the triple well could correspondingly 
deviate from the ideal sine-square form which we are using here. However our results and 
findings are to a large extent independent of the exact form the triple well potential and 
the boundary conditions, and similar effects would be detected for various experimentally 
realizable setups. 

The last term of H represents the ID contact-interaction potential. Under cylindrical 



harmonic confinement, the effective ID coupling becomes giD = ( 1 — ^~^^)~^ jl6| . 
where Oq is the free-space s-wave scattering length, and a± = (h/Muj±Y^'^ denotes the 
transverse-confinement length in terms of the frequency of the transverse harmonic confine- 
ment, u±. It is clear from the expression for Qid that the interaction strength can be tuned 
by either changing the scattering length oq, or the transverse confinement frequency uj±, 



appealing to field-induced Feshbach resonances 



13l-ll5|. or confinement induced resonance 




Figure 1. (a) Sketch of the triple weh potential as well as the first three Wannier levels in each well, 
(b) Response of the on-site energy of different number states to the interaction strength. The red 
line gives the response of the triple modes, and the blue and green lines provide the response of the 
single modes and pair modes, respectively, number states belonging to different bands are denoted 
by different line types: solid - ground band, dashed - first excited band, dash-dotted - second excited 
band, dotted - third excited band. Arrows indicate the resonance points of the triple mode with 
various modes. 



171, ll8|, respectively. 



In the following we rescale (1) in units of the recoil energy Er = h^K'^/2M, while the 
space and time are given in the unit of and h/Ept, respectively. This amounts to setting 
h = M = K, = 1. The coupling of the interaction potential now reads g = 2gi£)M/h'^K,. 
The triple well potential turns to Vtr = Vosin'^x, with hard wall boundary conditions at 
X = ±37r/2. In this work, Vq is large enough such that each well confines three localized 
single-particle Wannier states: the ground, first excited and second excited Wannier states. 
As a result of the hard wall boundary conditions, the on-site energy of the Wannier states 
in the left and right well e^^^j^ is slightly larger than that of the Wannier states of the same 
energetic level in the middle well ej^^ {i G [0, 2]), as shown in figure 1(a). 



B. Computational Approach: The Multi- Configuration Time-Dependent Hartree 
Method 



We apply the numerically exact Multi-Configuration Time-Dependent Hartree (MCTDH) 



method 
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j|-|5l|. a wave packet dynamical approach to the a&-initio solution of multi- 
dimensional time-dependent Schrodinger problems. In MCTDH the many-body wave func- 
tion is expanded in terms of Hartree products of single particle functions of a corresponding 
basis: 

jl,---,jN 

where ipj-{xi) are the single particle functions for the degree of freedom Xj. Substituting this 
Ansatz into the Schrodinger equation, ihdt^ = , leads to a set of differential equations 
for y4j^ jj^(t) and ipj^{xi,t). Integrating these differential equations, we obtain the time 
evolutions for arbitrary initial conditions \I'(xi, Xat; 0). MCTDH can also be used to 
calculate stationary states of given system by relaxation of a proper wave function, i.e. 
imaginary time propagation. The advantage of MCTDH compared to exact diagonalization 
is that it optimizes the single particle function ipj-{xi,t) at each time step, which leads to 
a reduction of the total number of necessary single particle functions in order to achieve 
convergence. When dealing with bosonic system, MCTDH can be simply modified to fulfill 
the permutation symmetry preserved by bosons. This can be done by symmetrizing the 
coefficients Aj, and using a single set of {(fj{x,t)} for all the degrees of freedom. 

We aim at the tunneling dynamics of initial states where all bosons are localized in the 
same well. Such an initial condition can be experimentally prepared e.g. by adiabatically 
applying an external potential which rearranges the energies of minimum of the triple well, 
and all bosons are then confined to the lowest well, for instance, a linear tilt potential can 
confine all the bosons to the left well, and a harmonic potential originating at the middle 
well will localize all the bosons to the middle well; consequently this external potential is 
removed instantaneously. To prepare the initial state for our numerical calculations, we 
apply artificial hard wall boundary conditions on each side of the corresponding well, and 
let the bosonic wave function relax to the ground state inside the well via imaginary time 
propagation in the framework of MCTDH. Subsequently we turn off the artificial hard walls, 
and explore the time evolution "^{xi, ...,xi\f]t) in the triple well potential via MCTDH. 
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III. GENERALIZED NUMBER-STATE REPRESENTATION AND SPATIO 
-ENERGETIC CHARACTERISTICS OF THE TRIPLE WELL 



Here we describe our number-state representations including the contribution of different 
bands. Tfiis yields an illustrative picture which helps us to analyze and interpret our nu- 
merical results obtained by MCTDH, and also serves as an explanatory tool. Especially the 
interband processes can be interpreted in the framework of this representation. We consider 
an ensemble of repulsively interacting bosons confined in the triple well potential as intro- 
duced in Sec. II. The triple-well potential is deep enough such that the Wannier functions 
of the bosons belonging to different wells have very small overlap for not too high energetic 
excitation. The wave function of the bosons can be expanded in a series of generalized 
number states which denote the distribution of all the bosons among the three individual 
wells: 

\'^)-J2^N\NL,NM,NR)i. (3) 

N,i 

The summation runs over all possible configurations N = {Nl, Nm, Nji) and the energetic 
excitations denoted by i. Nl, Nm and Nr are the number of bosons localized in the left, 
middle and right well, respectively, with N^ + Nm + Nr — N. We define the spatial 
representation of the number states as: 

(x|7V^, Nm, NR)i = 50V(x - r^v), (4) 

where S refers to the symmetrization operation among all the bosons, and x = (xi, ...,xn), 
= f^i'^ ® ^M*^ ® {{rL,fM,i^R} label the minimum positions of the three wells). The 
wave function 05v satisfies: 

(E -^^i + E ^(^^■) + - ew</>w(x), (5) 

i=i j=i 

where W{x) — Vosin'^{x) with hard wall boundary conditions at x = ±7r/2, is the local 
confining potential, and e^v is the ith on-site energy, with i labeling the level of energetic 
excitation. The interaction potential in (5) takes the form: 

= |[ E ^'^^ji ~ ^j^) + E ^^^j^ ~ ^^2) + E ^^^^1 - ^32)]- 

h7^j2&[l,NL] ji7^j2&[NL+l,NL+NM] ji9^j2&[N-Nn+l,N] 

(6) 
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A few remarks are in order here. Equation (5) is the constituting equation for the three 
subsets of bosons with particle number Nl, Nm and Nr- Within each subset all interactions 
are taken into account whereas no interactions are taken into account for bosons belonging 
to different subsets, which is because we divide the triple well for our definition of number 
states into three adjacent single wells with hard wall boundary conditions. That is why 
Vi can be divided into three terms according to (6). The number states {'x\Nl, Nm, Nii)i 
(see (4)) can be obtained from the wave functions 0^(x) of (5), which are located in the 
interval [— 7r/2, vr/2], by performing the appropriate shift of arguments contained in S(f)^j^f{'x— 
Tn), which leads to A''^;^, and Nji bosons localized in the left ([— 37r/2, — 7r/2]), middle 
([— 7r/2, 7r/2]) and right well {[tt/2,3tt/2]), respectively. 

Let us discuss the response of the energies e*^ to the interaction strength. As the inter- 
action strength rises from zero to infinity, the wave functions of the number states change 
from the form of non-interacting bosons to that of fermionized bosons 2l|, |22|, |32|, in which 



bosons with strong repulsive interaction in one dimension avoid spatial overlap and acquire 
fermionic properties. Correspondingly, the on-site energy increases from the sum of different 
single-boson Wannier energies and saturates to the corresponding value for non-interacting 
fermions. The latter can be described as = Cq + e/(A'i, Nm, Nji), where Cq is the value 
for non-interacting bosons, which is independent of the distribution of the bosons among 
the three wells, and increases with i. ef{NL, Nm, A"^) is the increase of the on-site energy in 
the fermionization regime compared to the non-interacting regime, which is dependent on 
the distribution of {N^, Nm, Nr), i.e., the more the bosons are localized to the same well, 
the higher e/(A^L, Nm, Nr) is. 

When the on-site energy difference of two number states becomes much smaller than the 
effective coupling between them, significant tunneling will take place between them. For 
instance, the on-site energy of number states with lower i, but higher degree of localization 
of the bosons can cross that of other number states at a particular interaction strength, 
resulting in an enhanced tunneling between the corresponding number states. As such 
tunneling directly results from the contributions of excited number states, from higher bands 
of the system, we refer to this enhanced tunneling as "interband" tunneling. 

Let us briefly provide a comparison of the number states introduced here and in other 
works focusing on higher band effects. In previous works, the wave function is expanded 
with respect to the number states of non-interacting bosons, and the interaction potential 
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is treated as a perturbation to these states. This means that these number states are just 
products of non-interacting single-particle Wannier functions. The perturbative treatment 
of the interaction potential limits such an expansion to the low-interaction regime. In our 
case the generation of number states intrinsically treats the interaction and the confinement 
potential, i.e., Vi{x.) and W{x) in (5), in an inseparable manner, so that it becomes valid in 
the complete interaction regime, and uncovers the interplay between the interaction potential 
and the confinement potential. For instance, the interband tunneling discussed here can be 
easily understood in the number-state basis developed here but not in the number-state 
basis of non-interacting bosons applied in previous works. 



IV. INTERBAND TUNNELING: OBSERVATION AND ANALYSIS OF MECHA- 
NISM 

Interband tunneling is a general phenomenon insensitive to the total number of bosons, 
given all the bosons are well localized in the individual wells, and the minimal system to 
unravel this phenomenon is that of three repulsively interacting bosons in the triple well. In 
the following we demonstrate this process at hand of the system of three bosons in a triple 
well. 

Firstly we briefiy discuss the number-state properties of the three-boson system, especially 
the properties of the on-site energies of the number states. The number states can be divided 
into three categories: the "single mode" {|1, 1, !)«}, implying that all the bosons are localized 
in different wells, the "pair mode" {|2, 1, 0)^, |2, 0, 1)^, |1, 2, 0)^, |0, 2, 1)„ |1, 0, 2)„ |0, 1, 2)^}, 
refering to that two bosons are localized in the same well, and the "triple mode" {|3, 0, 0)j, 
|0, 3,0)i, |0, 0,3)i}, which refers to all the bosons residing in the same well. The three 
categories of number states show a different response to the increase of the interaction 
strength, as schematically shown in figure 1(b). For the single modes, the on-site energy is 
barely affected by the interaction strength and remains constant in the complete interaction 
regime as the summation of the on-site energy of single-boson Wannier states in each well, 
i.e., e*^ = + + where the configuration {l,m,n) depends on the degree of 
energetic excitation i. The on-site energy of the pair modes always coincides with that of 
the corresponding single modes in the non-interacting regime, and saturates to the next 
upper band of the single modes for very strong interaction. For the triple modes, the on-site 
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energy saturates to the third upper band of the single modes, and crosses several energy 
levels of single and pair modes in the intermediate regime. The on-site energy spectrum of 
the number states is shown in figure 1 (b) , in which we can see the "band-like" structure of the 
number states. In the figure, the "energy gap" between "bands" is to a good approximation 
the on-site energy difference of Wannier states of adjacent energy levels, while the energy 
splitting within a "band" results from the on-site energy difference of the Wannier states 
in the outer and middle well, i.e., s'-j^jj^ — e^^. When the on-site energy of the triple modes 
crosses that of the single and pair modes of different energy levels, the interband tunneling 
process will take place, also marked in figure 1(b) with corresponding arrows. It is worth 
mentioning here that in the single band approximation, only number states with i = are 
considered, and the interband tunneling cannot be explained within this approximation. 

In the following we provide several examples of the interband tunneling process by nu- 
merically exact quantum dynamical studies. 

A. Single-boson tunneling 

In this section, we demonstrate the single-boson interband tunneling of three bosons 
initially locahzed in the same well of a triple well. Let's concentrate first on the case where 
all atoms are loaded in the left (or equivalently rightmost) well, i.e., the initial state is of 
the form |3,0, 0)o. The single-band model predicts self-trapping of all the bosons in the 
initial well when the interaction strength is large enough. Figure 2 shows the oscillation 
of the population in each well for the tunneling of three bosons initially localized in the 
left well for different interaction strengths. As shown in figure 2(a), for the interaction 
strength g — all the bosons remain essentially trapped in the left well, and tunnehng is 
strongly suppressed. This indicates that the bosons become self-trapped due to the repulsive 
interaction. In fact, the local Lieb-Liniger parameter 7 = g/Np{x) for the initial state can 
be estimated as 7 0.03 <C 1, signalling weakly interacting bosons, and an effective lowest- 
band model is applicable, with the local interaction energy Uq ~ 5.6 and a hopping parameter 
J ~ 0.01. Then we obtain Uq/J^ 560 1, which confirms the self-trapping phenomenon 
as shown in figure 2(a). 

In figure 2(b), however, for g — 3.26, where the single-band model predicts even stronger 
self-trapping than that oi g — 0.1, we observe the enhanced tunneling. At this value of 
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Figure 2. The population oscillations of three bosons in the triple well system, which are initially 
localized in the left well. The interaction strength is (a) g = 0.1 and (b) g = 3.26. Both figures 
demonstrate the emergence of enhanced tunneling windows on top of the suppressed-tunneling 
background, see discussion in text. 



g, the interaction parameter is 7 1.2, implying that on-site correlations are indeed non- 
negligible, and the single-band model is inapplicable here. In figure 2(b), the population in 
the left well oscillates between three and two, and that of the middle and right well oscillates 
between zero and one approximately. This implies that a single boson tunnels to the middle 
and right well with the other two bosons remaining localized in the left well. 

To clarify the mechanism underlying the enhanced tunneling in figure 2(b), we monitor 
the contributin from different number states, \i{NL, Nm, Nfi\'${t))\^, as a function of time. 
Here the number states [A^^^:, Nm, Nfi)i are computed numerically utilizing improve relaxation 
within MCTDH, and then |\I/(t)) is projected onto them at different times. The result is 
shown in Figure 3. Only the states |3,0,0)o, i.e. the initial state, and |2, 1, 0)i,|2, 0, l)i 
possess a significant contribution in the course of the tunneling process, indicating that the 
tunneling mainly takes place between these states. |2, l,0)i and |2,0,l)i refer to states 
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Figure 3. Probability of number states for the tunneling process shown in figure 2(b), for the states 
|3, 0, 0)0, [2, 1, 0)1 and |2, 0, among which the tunneling mainly takes place. For this value of g, 
these three number states are in resonance, as marked by the blue arrow in figure 1(b). 

where two bosons localize in the two-boson ground state of the left well while one is in the 
first excited Wannier state of the middle and right well, respectively. In this way a boson 
hops to the first excited Wannier state of the middle and right well during the tunneling 
process, and the tunneling is indeed a single-boson tunneling. 

Besides the number states occupation, the tunneling process can also be demonstrated 
by the one-body density evolution. Figure 4(a) shows the one-body density at different 
time instants of the tunneling process. The density profiles can offer additional evidences 
of the excitation into different bands throughout the time evolution. At t = the density 
profile is roughly a Gaussian localized in the left well, and this corresponds to the initial 
state |3, 0, 0)o. As time evolves, the density profile populates the middle and right well, and 
instead of the typical single-particle ground-state profile of a Gaussian, the profile in the 
middle and right well shows a nodal structure, which is typical for a single boson in the first 
excited Wannier state, thereby verifying the contribution of the states |2, 1, 0)i and |2, 0, l)i 
to which we referred above. The density profile appearing in the middle and right well is a 
spatial sign of the occupation of the first excited Wannier states in each of these wells, and 
based on this, we can give a schematic picture of this tunneling process as in figure 4(b), in 
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Figure 4. (a) The one-body density at various time instants. Initially all bosons are localized in 
the left well (black line), then one boson tunnels to the middle well (red line), and to the right 
well (blue line). The nodal structure in the middle and right well demonstrates that the boson in 
the middle and right well are in the first excited Wannier level, (b) Visualization of the tunneling 
process: in the tunneling process a boson tunnels forth and back between the left, middle and right 
well, with the other two bosons remaining in the left well. The dashed lines in the middle and right 
well illustrate the Wannier energy levels as in figure 1(a). The parameters are the same as in figure 
2(b) and figure 3. 

which a single boson initially localized in the left well tunnels among all the three wells, and 
the other two bosons remain localized in the left well. When the boson tunnels to the middle 
and right well, it actually occupies the first-excited Wannier level in the corresponding well. 

Now we proceed to demonstrate an alternative single-boson tunneling, in which the tun- 
neling boson can hop to even higher Wannier states at particular interaction strengths. As 
an example, we choose the initial state where all bosons are trapped in the middle well, 
i.e., |0,3,0)o, and the interaction strength g = 9.85 corresponding to 7 ^ 5.9. We again 
calculate the time evolution of the population in each well and the probability of different 
number states, as shown in Figure 5. This time the population in the middle well oscillates 
approximately in the interval [2, 3], while the population in the left and right well oscillates 
approximately within [0,0.5], as shown in figure 5(a). This indicates that a boson in the 
middle well tunnels to left and right well with the same probability. Figure 5(b) confirms 
this tunneling behavior with the probability evolution of the number states, in which the 
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Figure 5. (a) The population oscillation of three bosons in the triple well, for the interaction strength 
g = 9.85. The three bosons are initially localized in the middle well, (b) Time-dependence of the 
occupation of |0, 3,0)o, |1,2,0)3 and [0,2,1)3. As the system tunnels to the cat state |1,2,0)3 + 
10,2,1)3, one boson tunnels to the left and right well with the same probability, which gives rise 
to the coincidence of the population in the left and right well and that of the two number states 
11,2,0)3,10,2,1)3, as shown in (a) and (b), respectively. 

states |0, 3, 0)o, |1, 2, 0)3 and |0, 2, 1)3 share the major contribution and consequently tunnel- 
ing mainly takes place among these number states. |1, 2, 0)3 and |0, 2, 1)3 refer to the states 
of two bosons in the middle well with the third one in the second excited Wannier state of 
the left and right well respectively, and this indicates that the bosonic tunneling to the left 
and right well hops to the second excited Wannier state. 

To confirm and further analyze the above statement we show the corresponding one-body 
density for different times in figure 6(a). In the one-body density results, we see initially 
all bosons are trapped in the middle well. The deviation of this profile from a Gaussian 
indicates the onset of the fermionization process of the three strongly repulsively interacting 
bosons, which is similar to that of three strongly repulsively interacting bosons confined in 



the harmonic trap 
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53] . At later times, where a single boson tunnels to the left and 
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Figure 6. (a) The one-body density at various time instants for the tunnehng process in figure 
5. All bosons are initially localized in the middle well (black line), then one boson tunnels to the 
left and right well with the same probability (red and blue lines). The nodal structure in the left 
and right well indicates the occupation of the second excited Wannier states in both wells . (b) 
Visualization of the tunneling process: in the tunneling process a boson tunnels from the middle 
well to the outer wells. The dashed lines in the left and right well illustrate the Wannier energy 
levels. 

right well, there are two nodes in these wells demonstrating the occupation of the second 
excited Wannier energy levels in both wells. We also illustrate this tunneling process in figure 
6(b), and in this case the ensemble of bosons tunnels between the initial state |0, 3, 0)o, and 
the cat state |1, 2, 0)3 + |0, 2, 1)3. 

In conclusion, we demonstrated that tunneling emerging from a state of three bosons 
localized in the same well can be enhanced by the resonant coupling to number states 
relating to higher bands, and this resonant coupling gives rise to interband tunneling, in 
which a single boson tunnels to different excited Wannier states of different wells. 

The above-mentioned interband tunneling can also interplay with the external confine- 
ment to achieve a tunable tunneling not only to a certain band but also to a certain well. 
A slight tilt potential Vuh = 0.1 ■ x can e.g. be applied to the triple well, which detunes 
|2, 1,0)1, |2, 0,1)1, and now |3,0,0)o gets into resonance with |2,l,0)i and |2,0,l)i sepa- 
rately at different interaction strengths. Consequently a boson can selectively tunnel to the 
middle or right well at different interaction strengths. For instance, |3,0,0)o gets into reso- 
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Figure 7. Applying a tilt potential Vtj/f(x) = 0.1 ■ x, the tunneling of three bosons is shown with 
the initial condition that all bosons reside in the left well, for g = 3.0. The population in the left, 
middle and right well is provided. One boson exclusively tunnels to the middle well. 



nance with |2, 1, 0)i at the particular interaction strength g = 3.0, and consequently a single 
boson tunnels only to the middle well. Figure 7 provides the results of a tunneling process 
with tilt potential Vtm, in which the boson only tunnels to the middle well at g = 3.0. 



In experiments, the main difference of various realizations of the triple well potential is 
represented by the shift of the on-site energy of each well. The interband tunneling under 
tilt potential demonstrates that such tunneling is robust against the on-site energy shift, 
and consequently robust against deformations of the triple well potential from the ideal 
sine-square form. Earlier works also suggest that the tilt potential can be used to transport 
bosons in multiwell both dynamically [54,] and adiabatically j38|, and we point out that the 
interplay between the interaction potential and tilt potential can realize a specific transport 
of bosons not only into particular wells but also to particular Wannier states in corresponding 
wells. 
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Figure 8. Tunneling of three bosons initially localized in the middle well for g = 5.8. (a) Population 
oscillation of the left well (red line) , the middle well (black line) , and the right well (blue line) . The 
plots of the population in the left and right well practically lies on top of each other, indicating the 
correlated property of the two bosons tunneling to these wells, (b) Time occupation probability of 
the states |0,3,0)o, |1, 1, 1)6, |l,2,0)i and [0,2, l)i. 

B. Two-boson correlated tunneling and beyond 



An interesting phenomenon in case of a double well is pair tunneling 30|, l3l|, for which 
two repulsively interacting bosons tunnel as a pair forth and back. Pair tunneling arises since 
number states corresponding to bosonic pairs in each well are resonant with each other, and 
detuned from all other states. Pair tunneling is a special case of correlated tunneling, in 
which particles do not tunnel independently. When higher bands come into play, various 
types of correlated tunneling occur and can be addressed by simply tuning the on-site energy 
of the initial state into resonance with the desired state. 

As an example, we demonstrate numerically a correlated tunneling process, in which two 
of the three bosons initially localized in the middle well simultaneously tunnel to the left 
and right well, respectively. We take as the initial state |0, 3, 0)o for g = 5.8, where |0, 3, 0)o 
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Figure 9. (a) One-body density at different times for tlie tunneling as in Figure 8. Initially all the 
three bosons are localized in the middle well (black line) , and gradually two of them tunnel to the 
left and right well (red and blue lines), (b) illustration of such tunneling process. 

is resonant with |1, 1, 1)0, which refers to the configuration of a single boson located at the 
ground Wannier state of the middle well while the other two are in the first excited Wannier 
state of the left and right well, respectively. We show the evolution of the population in each 
well and the probability of number states in Figure 8. From the population oscillation shown 
in figure 8(a), we see that the occupation of the middle well oscillates within the interval 
[1,3], and each occupation of the left and right well oscillates within [0,1]. This indicates 
that two bosons initially localized in the middle well propagate in a correlated manner to 
the left and right well in the course of the tunneling process. From the time evolution of 
the probability of number states, shown in figure 8(b), we find that the tunneling mainly 
takes place between |0,3,0)o and |1, l,l)e, with a minor intermediate contribution from 
|1, 2, 0)1 and |0, 2, l)i. According to the configuration |1, 1, 1)6, this means that during the 
tunneling two of the three bosons initially localized in the middle well tunnel to the first 
excited Wannier states of the left and right well simultaneously. 

In figure 9(a) we show the corresponding one-body density for different times to spatially 
resolve the interband tunneling process. The nodal pattern of the density profile in the left 
and right well confirms the occupation of the first excited Wannier states in these wells. 
Figure 9(b) is a schematic plot of such a tunneling process, based on the discussion above. 

Correlated tunneling refers to the notion of tunneling of interacting particles that don't 
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move independently, in comparison with the independent tunnehng of free particles. In the 
single band approximation it is impossible to tune arbitrary number states into resonance via 
a change of the interaction strength, which is due to the different on-site energy dependence 
of number states on the interaction strength. This, however, becomes possible when we 
take into account higher bands, and in principle tunneling between arbitrary spatial and 
energetic configurations is achievable. As an example, we have shown here the correlated 
tunneling of two bosons initially in the middle well simultaneously to the left and right well. 

Such interband tunneling can be straightforwardly extended to larger systems containing 
more bosons. We take four bosons in the triple well for instance. In such a system, there 
are four different categories of number states. In addition to the single mode, pair mode and 
triple mode, there exists the "quadruple modes" of {|4, 0, 0)j, |0, 4, 0)j, |0, 0, 4)j}, referring to 
four bosons in the same well. We can further separate the pair mode into two categories: the 
"double-pair modes", such as |2, 2, 0)j, referring to the case that the four bosons are divided 
into two pairs and each of the pairs occupies a different well, and the "single-pair mode", 
containing a pair and two separate bosons, such as |1, 2, Now in this system, besides the 
single-boson tunneling and two-boson correlated tunneling in the three-boson case, we can 
even realize correlated tunneling of pairs, by tuning the number states |0, 4, 0)o and |2, 0, 2)i 
into resonance, for instance. 



V. SUMMARY AND CONCLUSIONS 

We demonstrate that in a system consisting of a small ensemble of bosons confined in a 
one-dimensional triple well, several windows of enhanced tunneling are opened in the strong 
interaction regime, where in general the background of suppressed tunneling dominates. 
Such enhanced tunneling results from resonant coupling between the initial state and various 
excited number states, and can be only understood within an exact treatment of the system 
since it involves many bands, which is why we name it interband tunneling. When the initial 
state resonantly couples to different excited number states, various tunneling processes can 
be realized, which manifests itself as a potential tool for controllable dynamical transport 
of bosons. As an example, we demonstrate the single-boson tunneling to the first and 
second excited Wannier states, and the two-boson correlated tunneling in the system of three 
repulsively interacting bosons in the triple well with numerically exact quantum dynamical 
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studies. The interband tunneling we discuss here can be straightforwardly generalized to 
multiwell systems with more bosons. In this way, tunneling between arbitrary spatial and 
energetic configurations of bosons in a multiwell trap are in principle achievable just by 
tuning the interaction strength, and as a consequence controllable dynamical transport of 
bosons to specific wells and specific in-well energy levels becomes possible. This also opens 
the doorway to an interpretation and controllable preparation of the general non-equilibrium 
quantum dynamics in multiwell systems. 

During the interband tunneling process, bosons tunnel to excited number states, and 
this gives the opportunity to couple the bosons to photons via induced emission processes. 
Such coupling can map the properties of number states, such as the spatial configuration, 
to emitted photons, and may manifest itself as a possible controllable photon source. 

To understand the interband tunneling, we introduce a basis of generalized number states, 
which treat the interaction in a non-perturbative way such that these are valid even in the 
strong-interaction regime. Our number states are only meaningful and defined for large Vq. 
A main drawback of our method is that the interacting number states cannot be analyti- 
cally generated, and we have to rely on a numerical approach to obtain them. However this 
cannot be considered as a serious drawback, since the gain in terms of analysis and trans- 
parent interpretation of the dynamical process is substantial, allowing us to understand and 
consequently design the non-equilibrium tunneling dynamics in optical lattices. We also 
notice that efforts have been put forward to provide an analytical description in the strong 



interaction regime 
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